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1 Introduction 

The production of energetic photon pairs in hadronic collisions is a valuable testing ground 
of perturbative quantum chromodynamics (pQCD). The emission of a pair of photons from 
hard parton-parton scattering constitutes a particularly clean test of perturbation theory in the 
collinear factorisation fV, ^ and kj factorisation [3J approaches, as well as soft-gluon logarith- 
mic resummation techniques [4J. A comprehensive understanding of photon pair production 
is also important as it represents a major background in certain searches for rare or exotic pro- 
cesses, such as the production of a light Higgs boson, extra-dimension gravitons, and some 
supersymmetric states. 

This paper presents a measurement of the production cross section for isolated photon pairs in 
proton-proton collisions at a centre-of-mass energy of 7 TeV, using the Compact Muon Solenoid 
(CMS) detector at the Large Hadron Collider (LHC). Isolated photons produced in the hard 
scattering of quarks and gluons are henceforth referred to as signal photons and the remaining 
photons as background photons. A pair of signal photons will be referred to as a diphoton. The 
data sample was collected in 2010 and corresponds to an integrated luminosity of 36.0 pb^^. Re- 
cent diphoton cross-section measurements have been performed by the DO |3 and CDF [61 |Zl 
Collaborations in proton-antiproton collisions at a/s = 1.96 TeV, and by the ATLAS Collabora- 
tion at the LHC |i8J. 

The CMS detector consists of a silicon pixel and strip tracker surrounded by a crystal electro- 
magnetic calorimeter (ECAL) and a brass/ scintillator sampling hadron calorimeter (HCAL), all 
in an axial 3.8 T magnetic field provided by a superconducting solenoid of 6 m internal diam- 
eter. The muon system is composed of gas-ionization detectors embedded in the steel return 
yoke of the magnet. In addition to the barrel and endcap detectors, CMS has an extensive 
forward calorimetry system. A more detailed description of CMS can be found elsewhere [9J . 

In the CMS coordinate system, 9 and (p respectively designate the polar angle with respect to the 
counterclockwise beam direction, and the azimuthal angle, expressed in radians throughout 
this paper. The pseudorapidity is defined as f/ = — In [tan |] . 

Distance in the (?/, (p) plane is defined as R = ^y {Arj)^ + (^<p)^- The transverse energy Ej of 
a particle is defined as Ej = E sin 9, where E is the energy of the particle, and the transverse 
momentum is pT = psm9. The rapidity is defined as y = jhn [§z^]/ with pz being the 
longitudinal momentum with respect to the beam axis. 

The electromagnetic calorimeter, which plays a major role in this measurement, consists of 
nearly 76 000 lead tungstate crystals. It is divided into a central part (barrel) covering the re- 
gion < 1.48 and forward parts (endcaps) extending the coverage up to |?/| < 3 for a particle 
originating from the nominal interaction point. The crystals are arranged in a projective ge- 
ometry with a granularity of 0.0174 in both the rj and cp directions in the barrel, and increasing 
with r] from 0.021 to 0.050 in the endcaps. A preshower detector, consisting of two planes of 
silicon sensors interleaved with 3 radiation lengths of lead, is placed in front of the endcaps to 
cover the pseudorapidity region 1.65 < \}]\ < 2.6. 

The differential cross section is measured as a function of variables that are particularly relevant 
in searches for rare processes or to characterise QCD interactions (e.g. O): 

• the diphoton invariant mass, m^y; 

• the azimuthal angle between the two photons, Acpy^; 

• the photon pair transverse momentum, pT,^^ = J Pt,7i^ + Pt,72^ + 2 PT,yiPTy2 Acp^y, 
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where pT,yi and pT,y2 the magnitudes of the transverse momenta of the two pho- 
tons; 

• I cos 9* I = |tanh \, with Ay^y being the difference between the two photon ra- 
pidities. At lowest order in QCD, 6* is the center-of-mass scattering angle for the 
qq 77 and gg 77 processes. 

The event selection requires at least one isolated photon with Ej > 23 GeV and a second iso- 
lated photon with Ej > 20 GeV, separated by R > 0.45. The measurements are performed in 
two pseudorapidity regions, one with < 1.44, and the other defined by the tracker accep- 
tance \tj\ < 2.5, but excluding the transition region between the barrel and endcap calorimeters, 
1.44 < 1?^ I < 1.57. For convenience the widest ?/ range without the transition will be referred to 
as |f/| < 2.5 throughout the paper. 

The asymmetric thresholds on the photon transverse momenta avoid the infrared sensitivity 
affecting the fixed-order calculations liTOllTTI and simplify the comparison of the measurements 
with the theoretical predictions. 

All simulation results are based on the PYTHIA 6.4.22 fT2] event generator, with the z2 tune, the 
CTEQ6L parton distribution functions (PDFs) LISJ, and a GEANT4 [14] modelling of the detector. 
The z2 tune is identical to the zl tune described in ll5| except that z2 uses the CTEQ6l PDFs 
while zl uses CTEQ5l ||T6| . At the generator level, a prompt photon is considered as signal if 
the sum of the generated transverse momenta of all the particles within a cone R < 0.4 around 
the photon direction is less than 5 GeV. 

Event selection and background discrimination are presented in Sections |2] and |3] The determi- 
nation of the signal yield and the measurement of the cross section are described in Sections |4] 
and|5] Systematic uncertainties are detailed in Section|6] Results are discussed in Section|8]and 
compared with the theoretical predictions introduced in Section [7[ 

2 Event Selection 

Photon candidates are reconstructed by clustering the energy deposited in the ECAL (iTl [181 
crystals. CMS is equipped with a versatile trigger to adapt to the steady increase in the LHC 
instantaneous luminosity. In this measurement, three trigger settings were used for three suc- 
cessive data-taking periods. They require two photon candidates, with a threshold of either 
15 GeV or 17 GeV on the transverse energy. For the last period, with the highest instantaneous 
luminosity, a weak isolation requirement is applied on one of the two photon candidates. For 
the three periods, the trigger efficiency for events passing the analysis selections described in 
the following paragraphs is estimated from simulated events to be greater than 99.9%. The 
offline event selection requires one photon candidate with Ej > 23 GeV and a second photon 
candidate with Ej > 20 GeV, each within the fiducial region defined in the introduction. The 
candidates are required to be separated by R > 0.45 to avoid energy deposits related to one 
candidate overlapping with the isolation region of another candidate. 

Photon identification criteria requiring the deposits in the calorimeters to be consistent with an 
electromagnetic shower are applied to the two candidates. The criteria are based on the spread 
along rj of the energy clustered in the ECAL, henceforth referred to as CTrjtj, and on the ratio H/E 
of the energies measured in the HCAL and ECAL (loose selections of Ref. fflSij ). 

The photon candidates are required to be isolated. The sum of the transverse momenta of 
charged particles measured by the tracker and the sum of the transverse energy deposits in 
the HCAL, both defined within a cone of radius R = 0.4 around the photon direction, must 
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each be less than 2 GeV in the barrel and 4 GeV in the endcaps. HCAL deposits in a cone of 
radius R = 0.15 are excluded from the sum, as well as tracks in a cone of radius R = 0.04 
and within a strip of Arj = 0.03 along the cp direction, which can potentially contain tracks 
of an electron-positron pair from the conversion of a photon in the tracker material. The sum 
of the transverse energy deposited in the ECAL in a cone of radius R = 0.3 is required to 
be less than 20% of the photon transverse energy, in order to be consistent with the online 
trigger requirements. Excluded from the sum is the energy deposited within a cone whose 
radius corresponds to 3.5 crystals along rj and within a 5-crystal-wide strip extending along 
the cp direction. In addition, we require that no charged particle with the following properties 
impinge on the ECAL within a cone of radius R = 0.4: transverse momentum pj > 3 GeV, 
impact parameters with respect to the primary vertex in the transverse and longitudinal planes 
of less than 1 mm and 2 mm, respectively, and one associated hit in the innermost layer of the 
pixel detector. Tracks corresponding to such particles are henceforth called impinging tracks. 
The electron contamination is further reduced by imposing an additional veto on the presence 
of hits in the layers of the pixel detector along the direction of the photon candidate. 

3 Signal and Background Discrimination 

The photon candidates in the selected event sample are designated as signal photons, back- 
ground photons from hadron decays (most of which are misidentified pairs of collinear photons 
coming from neutral meson decays), or misidentified electrons. The background to diphoton 
pair events is thus made up of photon+jet and multijet events, with respectively one and two 
background photons from neutral hadron decays, and Drell-Yan events, with two misidenti- 
fied electrons. 

The contamination from Drell-Yan events is estimated from simulation using the next-to-leading- 
order (NLO) POWHEG generator ||T91-|2T|. which agrees well with our own Drell-Yan measure- 
ment II22I . The diphoton cross-section measurement is corrected for this contamination, which 
amounts to about 12% in the diphoton mass range 80-100 GeV around the Z peak. This proce- 
dure has a negligible impact on the systematic uncertainties. 

Background photons from photon+jet and multijet events are produced in jets alongside other 
particles, which tend to widen the deposits in the ECAL. An isolation variable X based on the 
energy in the ECAL is used to statistically estimate the fraction of diphoton events among the 
selected candidates. This variable is constructed to minimise the dependence on the energy 
deposited by minimum-ionising particles (MIPs) such that its distribution for the background 
can be obtained from the data by means of the impinging-track method described below. It is 
defined as the sum of the transverse energy of the ECAL deposits with Ej > 300 MeV (MIP 
veto), within a hollow cone centred on the photon impact point, with an inner radius of 3.5 
crystal widths and an outer radius of R = 0.4. Deposits assigned to the photon itself or falling 
within a 5-crystal-wide strip extending along cp and centred on the photon impact point are 
removed. Thus, deposits from photons converting into electron-positron pairs in the tracker 
material and spread along the cp direction do not contribute to the value of I. The variable I 
differs from the ECAL isolation used in the selection described in Section |2l 

As the distribution of I is different for signal photons and background photons, this variable 
can be used in a maximum-likelihood fit to extract the number of signal events in the entire 
selected sample. Figure [l] shows the probability density function of I, which was extracted 
from data with the methods described below. 

Contributions to the value of I for signal photons come from pileup (multiple proton-proton 
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collisions in the same bunch crossing) and underlying-event activity (multiple parton interac- 
tion and beam remnants from the same proton-proton collision). Since these contributions are 
independent of cp, the ECAL isolation probability density function /(X) is estimated from ran- 
dom cones using events with at least one isolated photon candidate. The value of X is calculated 
in a cone of radius R = 0.4 around an axis at the same value of as the photon candidate and 
at a random value of q) within a window of width n/l centred on the axis perpendicular to the 
photon direction, with the same exclusions applied to photon signals. The cone is required not 
to include any photon or electron candidates or jets. The function /(X) for signal photons is 
validated with two additional independent methods. Both methods exploit e^ and e^ from Z 
and W boson decays that do not radiate significantly in the tracker material. The e+ and e^ are 
selected with a constraint imposed on the fraction of bremsstrahlung energy emitted from the 
interaction in the tracker material. Such electrons and positrons leave ECAL energy deposits 
consistent with those of photons, and have a similar probability density function for I. The 
Z — )> e+e^ events are selected with stringent requirements on the identification criteria of the 
lepton pair and on its invariant mass, and the /(X) distribution is obtained directly from both 
leptons. In W — > eVg events, /(X) is obtained by exploiting the sPlot technique II23II . The miss- 
ing transverse energy projected along the lepton axis is used to estimate the probability of an 
event to be signal (W — ?► eVg) or background (Z — > e+e^, W — ?► rvj, 7 + jet(s), and QCD multijet 
processes). The value of X for the selected candidates is weighted accordingly to estimate the 
distribution of X. The uncertainty on /(X) is taken as the maximum difference between the 
distributions extracted from random cones and from electrons in Z and W events. In simulated 
events, the difference between /(X) for signal photons and for random cones is smaller than 
the uncertainty determined from data. 

For background photons, /(X) is extracted from a sample with less than 0.1% of signal-photon 
contamination. The sample is obtained by selecting photon candidates with one and only one 
impinging track. A cone of radius R = 0.05 around the track is excluded from the isolation 
area to avoid counting the energy deposited by the charged particle. The isolation variable 
X is then rescaled to take into account this additional exclusion. To validate this method, the 
X distribution is also extracted from a sample of events with two impinging tracks, one of 
the two being excluded in the computation of X. The latter distribution is compared to that 
obtained with the one-impinging-track sample, using the normal definition of X, i.e., including 
the energy deposits in the vicinity of the track. The agreement is within one standard deviation 
for the entire range of the X distribution, and the difference is taken as a systematic uncertainty 
on /(X) for background photons. 

The distributions /(X) show a moderate dependence on rj and on the pileup conditions, the 
latter being quantified by the number Wytx of primary vertices in the events (2.4 on average). 
The background distribution /(X) also depends on the transverse energy Ej of the candidate. 
Therefore, events in the sample used for the extraction of /(X) are weighted to reproduce the 
distributions of tj, n^tx, and Ej of the diphoton sample used for the cross section measurement. 
The effect of using the distributions from the diphoton sample to correct the biases in the back- 
ground and signal shapes used in the maximum-likelihood fit is addressed in the systematic 
uncertainty section. Section |6] 

4 Signal Yield Determination 

The number of diphoton events is obtained from a binned maximum-likelihood fit to the dis- 
tributions of the ECAL isolation variables of the two photons, Xi and X2, whose ordering is 
chosen randomly. Events are separated into three types: signal events (77) if both photons 
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are signal photons, background events with a signal photon and a background photon, and 
background events with two background photons. 

The likelihood function L that is maximised in the fit is 

N'°' N 3 

^=^nEw^m.4). (1) 

where N is the number of selected events, Nt is the number of events estimated in the fit for 
event type t, N*°* is their sum, and ji(Xx,X'i) is the probability for the ECAL isolation variables 
of the two photons to have values X\ and X2 for a given event type i. 

The probability density functions /f(X{,X|) for the three event types are obtained by multi- 
plying the probability density functions fiX) for single-photon candidates, assuming the two 
statistical variables X\ and X2 to be independent. Correlations between these two variables 
have been checked with simulation and are negligible. 

A total of 5977 events pass the selection criteria described in Section|2] These events are divided 
into three subsamples depending on whether both photons are in the barrel (2191 events), one 
is in the barrel and the other in the endcaps (2527 events), or both are in the endcaps (1259 
events). The fit is performed separately for each of the three subsamples and each bin of the 
four observables. An example of the fit for one bin of the wj-y-y spectrum is shown in Fig. |2]for 
events with both photons in the barrel ( | ?/ 1 < 1 .44) . 

The maximum-likelihood method is known to be biased for samples with small numbers of 
events [24J. This bias is estimated with Monte Carlo pseudo-experiments and the fit results are 
corrected for it. It is less than 10% of the statistical uncertainty for 80% of the fits and never 
exceeds half the statistical uncertainty. 

5 Cross-Section Measurement 

The differential diphoton cross-section measurement du I dX, for the variable X in the interval 
X„ is 

^(x) = (2) 

dX^ CAXiC{Xi)' ^ ' 

where N^ry is the number of signal events obtained from the fit, unfolded for the detector reso- 
lution and corrected for the Drell-Yan contamination; C is the integrated luminosity, AX, is the 
interval width; C is a correction factor for the effects of the detector resolution on the acceptance 
and on the efficiencies of photon reconstruction and identification. 

The number of signal events is unfolded || 25| for the detector resolution by inverting a response 
matrix T for each of the observables nijy, ^7,77/ A^-y-y, and |cos0*|, obtained from simulated 
events passing the selection requirements. The matrix elements T'^ are the probabilities of a 
selected event with the generated value of X in bin Xj- to be reconstructed with a value of X in 
bin X,. For a given interval Xj., the number of events after unfolding is related to the observed 

numbers of events in the different intervals X, by N}^^{Xk) = (T"^)'^'N-y-y(X,). Here, N-y-y(X,) 
is the signal yield corrected for the Drell-Yan contamination, as described in Section |3] Given 
the excellent energy resolution of the ECAL and the bin sizes, the matrix T is nearly diagonal, 
and thus no regularisation is applied in the unfolding procedure. 
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5 Cross-Section Measurement 
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Figure 1: Probability density functions of the ECAL isolation variable X for signal photons 
(solid blue) and background photons (dashed red) in the barrel (left) and in the endcap regions 
(right). 
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Figure 2: Fit to the photon ECAL isolation (Xi,X2) in the bin 100 < m^^ < 140 GeV for photons 
with \ tj\ < 1.44. The distribution of the isolation variable Xi of one photon candidate, arbitrarily 
chosen as the "first photon" and denoted with subscript "1", is displayed in the left figure, 
together with the fit result, integrated over I2; the shaded region shows the signal distribution, 
the dashed line represents the background contribution, while the solid line is the sum of the 
signal and background contributions. The same distributions for the second photon candidate 
are shown in the right figure. In this mass bin, the number of signal events is 72 ± 14, out of 
the total number of 161 selected candidates. 
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The correction factor C(X,) is defined as 



ATsim (Y) cdata 

nfY ^ recov '/ /-o-v 



where 



Nre^(X;) is the number of simulated events passing all the selection criteria, with generated 
values of X within the interval X,; 

Ngg^(X,) is the number of simulated events within the acceptance defined at the generator 
level (Section [T]), with generated values of X within the interval X,; 

gdata jg ^Yxe efficiency of the photon identification criteria measured from data; 

gsim jg ^Yxe efficiency of the photon identification criteria obtained from simulated events using 
the same technique as for e'^^'^. 



The efficiencies e*^^'^ and to observe a diphoton candidate are taken as the square of the 
efficiencies to observe a single photon. 

The efficiency for the requirements on isolation, cr^,j, and H/E is estimated with a "tag-and- 
probe" method l|26ll applied to a Z — > e+e^ sample selected from the full 2010 dataset. One 
lepton, the tag, is selected with tight reconstruction and identification criteria [27], while the 
other, the probe, is selected by requiring a constraint on the invariant mass of the lepton pair. 
The probes constitute a sample of unbiased electrons and positrons. The same constraint as 
discussed in Section |3] is applied on the fraction of bremsstrahlung energy emitted by the e+ 
and e^ interacting in the tracker material. This requirement ensures that the electromagnetic 
deposits of these electrons and positrons are consistent with those of a photon shower. The 
efficiency is computed by applying the requirements on isolation, cr^^, and H/ E to this sample, 
and then measuring the fraction of probes passing the selection. 

The efficiency for the requirement to have no impinging tracks within the isolation cone is 
estimated from data, from a control sample built using a random-cone technique on events with 
a single photon selected according to the identification criteria described above. The random- 
cone definition is that introduced in Section [s] for the extraction of /(X). Particles within the 
random cone hence come mainly from pileup and the underlying event. Quantities such as the 
number of impinging tracks or energy deposits in the isolation area are therefore assumed to 
have the same distributions as for isolated photons. The efficiency of the requirement to have 
no impinging track within the isolation cone is given by the ratio of the number of random 
cones passing this criterion to the total number of random cones. The efficiency of the veto on 
pixel hits is obtained from simulation. It is included in the N^eJ^/N|™ term of Eqn. (Sjl. 

The correction factor C is (80.8 it 1.9)% for the integrated cross section in the region \rj\ < 1.44, 
and (76.2 ± 3.3)% in the region \ri\ < 2.5. 



6 Systematic Uncertainties 

The systematic uncertainty on the reconstructed photon four-momenta is dominated by the 
ECAL energy scale, known to 0.6% in the barrel and 1.5% in the endcaps f2E\. The energy scale 
affects the value of the acceptance and induces bin-to-bin migrations in the differential cross 
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sections. The effect on the acceptance is relevant only in kinematic regions near the photon pj 
thresholds and results in an uncertainty of 40% in the most affected region, the lowest values 
of niyj. The uncertainty from the bin-to-bin migration is about 1%. 

The systematic uncertainty on the measured photon identification efficiency (e'*^'^ in Eqn. (jsj) 
is estimated by applying the tag-and-probe and random-cone methods on simulated events. 
The difference between the efficiency value obtained with these methods and that given by 
the fraction of simulated events passing the identification criteria is taken as the systematic 
uncertainty. The uncertainty from the acceptance and efficiency correction factor C is taken 
as the quadratic sum of the statistical uncertainties on the different factors of Eqn. (jsjl and the 
systematic uncertainty mentioned above. The systematic and statistical uncertainties on e*^^*^ 
total 1.9% for diphotons in the barrel and 3.3% for all diphotons. 

The systematic uncertainties on the signal and background isolation probability distributions 
/(X) are estimated with Monte Carlo pseudo-experiments in which /(X) is varied. The vari- 
ations correspond to the differences between the shapes of the nominal and validation distri- 
butions observed in the validation of the random-cone and impinging-track methods (Section 
|3](. In the first bin of the probability density functions, they are of the order of ±0.01 for the 
signal, and range from ±0.03 to ±0.05 for the background. The uncertainty on /(X) from its 
dependence on the distribution of photon transverse energy Ej, photon pseudorapidity t], and 
number of vertices nvtx is estimated from the change in /(X) when using the Ej, t], and nytx dis- 
tributions from the diphoton simulation instead of those from the diphoton event candidates 
in data. This contribution to the uncertainty is negligible. The overall systematic uncertainty 
from the /(X) distributions on the integrated cross section is about 8%, and varies from 4 to 
27% on the differential cross sections, depending on the bin and the subsample. 

A 4% uncertainty is assigned to the integrated luminosity ||29|. The various contributions to 
the systematic uncertainties are summarised in Table [l] 

Table 1: Contributions to the systematic uncertainties on the measured differential cross sec- 
tions for two pseudorapidity ranges. The uncertainties are computed for each bin of Figs. |3] 
to 10 The values listed below are averages. 



Uncertainty source 


\f]\ < 1.44 


<2.5 


Energy scale 
on acceptance 


1.5% 


2% 


Energy scale 

on bin-to-bin migration 


1% 


1.5% 


Signal and background 
distributions for /(X) 


7% 


9% 


Acceptance and efficiency 
correction factor C 


2% 


3% 


Luminosity 


4% 


4% 


Total 


8% 


11% 



7 Theoretical Predictions 



This section introduces the theoretical calculations whose predictions are compared to the ex- 
perimental data in Section |8] The leading contributions to the production of pairs of prompt 



9 



photons in pp collisions are the quark-antiquark annihilation (qq —> 77), gluon fusion (gg —> 
77), and gluon-(anti)quark scattering (qg — > 77q) processes. One or both photons come either 
directly from the hard process or from parton fragmentation, in which a cascade of successive 
collinear splittings yields a radiated photon. Contributions from the quark-antiquark annihila- 
tion process and the single- and double-fragmentation processes are calculated up to order oisoi^ 
with the DIPHOX 1.3.2 program [IJ. The contributions from the gluon fusion process, includ- 
ing the one-loop box diagram of order oi'^oc'^, the interference between the one- and two-loop 
box diagrams, and the real emission one-loop "pentagon" gg — > 77g, both of order u^a'^, are 
calculated with the GAMMA2MC 1.1.1 program 0. The fragmentation function BFG set II fSO} 
is used in the calculation. Although they are higher-order processes, the gluon fusion contribu- 
tions are quantitatively comparable to those from quark-antiquark annihilation in the diphoton 
mass range of interest, due to the significant gluon luminosity in this mass range at the LHC. 
The three theoretical scales, renormalisation, initial factorisation, and fragmentation, are set to 
the diphoton mass value. 

The photons are required to be within the kinematic acceptance defined in Section [l] An ad- 
ditional isolation requirement at the parton level is imposed by requiring the total hadronic 
transverse energy deposited in a cone of radius R = 0.4 centred on the photon direction to 
be less than 5 GeV. Particles resulting from underlying-event activity and hadronisation are 
not included in partonic event generators such as DiPHOX and GAMMA2MC. The fraction of 
diphotons not selected due to underlying hadronic activity falling inside the isolation cone is 
estimated using the pythia 6.4.22 [El event generator with tunes z2, d6t L3U, pO [32], and 
DWT [31J. A factor of 0.95 ± 0.04 is applied to the parton-level cross section to correct for this 
effect. 

The uncertainties associated with parton distribution functions and the strong coupling con- 
stant Ois are determined according to the PDF4lhc recommendations [33J. The diphoton cross 
section is computed with three different PDF sets (CTlO [34J, MSTW08 [35|, and NNPDF2.1 |36|), 
taking into account their associated uncertainties and the uncertainties on a,. The respective 
preferred as central value of each PDF set is used, and as is varied by ±0.012. The value for 
the cross section is taken as the midpoint of the envelope of the three results, including the 
uncertainties (68% confidence level envelope). The uncertainty on the cross section is taken to 
be the half-width of the envelope. 

The theoretical scale uncertainties are estimated by varying the renormalisation, initial factori- 
sation, and fragmentation scales by factors of 1 / 2 and 2, keeping the ratio between any two 
scales less than 2 (for example the combination 0.5 niyy, Inijy, nijy is not considered). The 
uncertainty is taken to be the maximum difference in the resulting cross sections. 

8 Results 

The integrated diphoton cross sections obtained for the acceptances defined in Section [l] are 
(7(pp ^ 77)||,,l<i.44 = 31.0 ±1.8(stat.) +|° (^yst-) ± 1.2 (lumi.) pb , 

c^iPP ^ 11)\m<2.50 = 62.4 ±3.6 (Stat.) (syst.) ±2.5 (lumi.) pb. 

The theoretical calculation described in the previous section predicts 

(r{pp ^ 17\\<iM = 27.3 13.0 (g^^j^g) ^ ^ ^ (p^p) ^ 

<^{PP ^ 17)\\r,\<2.50 = 52.7 ;5.8 (gcaies) ±2.0 (PDF) pb. 
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9 Summary 



The integrated cross sections obtained from the calculation are consistent with the measure- 
ments within the experimental and theoretical uncertainties. 

The differential cross-section measurements as functions of niyy, A<pjy, Pt,77/ and |cos0*| for 
the two pseudorapidity ranges are shown, along with the theoretical predictions, in Figs. |3] to 



10 The 4% uncertainty on the integrated luminosity is not included in the error bars. The val- 
ues of the cross sections are given in Tables|2]to|5] As can be seen in Figs.[7|and|8j the theoretical 
predictions underestimate the measured cross section for Acpy^ < 2.8. In the leading-order (LO) 
diagrams of gluon fusion and quark-antiquark annihilation 2 — > 2 processes, the two photons 
are back-to-back because of momentum conservation. Therefore, the LO term does not con- 
tribute to this phase space region, which thus only receives contributions from NLO terms for 
both the direct and fragmentation diphoton production processes. 

The contribution for Acp^^ < 2.8, combined with the requirements of Ej > 20 and 23GeV on 
the two photons, is responsible for the shoulder around 40 GeV in the diphoton differental pr 
distribution of Figs. |5]and|6] This contribution also populates the region below 30 GeV in the 
diphoton mass distribution shown in Figs. |3] and [4j In these two regions of the pT,yy and niyy 
spectra, the theoretical cross section is lower than the measurement, consistent with the deficit 
for Afjy < 2.8. 

Comparison of the measurements of the |cos 6* \ spectra with theoretical predictions, shown in 
Figs.|9]and 10 reveals an underestimation from the theory at large |cos 9* \ values, which is more 
significant for the central rapidity range < 1.44). Similar discrepancies have previously 
been observed in diphoton production at hadron colliders ||5H8H37| as discussed in Ref. II38II . 



Table 2: Measured diphoton differential cross section as a function of niyy for the two photon 
pseudorapidity ranges, with statistical (stat.) and systematic (syst.) uncertainties. 



da/druyj [pb/GeV] 





[GeV] 




\v\ 


< 1.44 






\v\ 


< 2.5 










Stat. 


syst. 




Stat. 


syst. 


0- 


-30 


0.0299 


±0.0071 


±0.0069 


-0.0086 


0.050 


±0.013 


±0.014 


-0.024 


30- 


-40 


0.061 


±0.030 


±0.015 


-0.018 


0.127 


±0.049 


±0.035 


-0.061 


40- 


-45 


0.097 


±0.088 


±0.020 


-0.020 


0.28 


±0.17 


±0.06 


-0.07 


45- 


-55 


0.77 


±0.12 


±0.06 


-0.05 


1.40 


±0.20 


±0.14 


-0.12 


55- 


-65 


0.70 


±0.10 


±0.05 


-0.04 


1.43 


±0.18 


±0.10 


-0.09 


65- 


-80 


0.408 


±0.059 


±0.030 


-0.031 


0.80 


±0.11 


±0.07 


-0.06 


80- 


-100 


0.175 


±0.031 


±0.013 


-0.012 


0.365 


±0.063 


±0.041 


-0.037 


100- 


-140 


0.070 


±0.012 


±0.003 


-0.003 


0.142 


±0.028 


±0.020 


-0.018 


140- 


-200 


0.0102 


±0.0035 


±0.0007 


-0.0006 


0.054 


±0.015 


±0.006 


-0.006 


200- 


-300 


0.0022 


±0.0011 


±0.0001 


-0.0001 


0.0084 


±0.0060 


±0.0023 


-0.0019 



9 Summary 

The integrated and differential production cross sections for isolated photon pairs have been 
measured in proton-proton collisions at a centre-of-mass energy of 7 TeV, using data collected 



11 



CMS \/s = 7 TeV, L = 36 pb"' 

I < 2.5, Ej^ > 20, 23 GeV 

— I I I I I I I I I I I I I 

Theory DIPHOX + GAMMA2MC- 
Measured 
Stat, uncertainties 
Stat. © syst. uncertainties 
Ttieoreticai scaie uncertainties 
PDF + a„ uncertainties 



CMS 



Vs = 7 TeV, L = 36 pb"' 




10"' 



CO 
T3 



4 

3.5 
3f 
2.5 

2r 
1.5r 
1 

0.5 




|ti| < 2.5, E-r^ > 20, 23 GeV 



-0.5 



(data - theory) / theory 
Stat, data uncertainties 
Stat. © syst. data uncertainties 
Theoretical scale uncertainties 
I PDF + uncertainties 



50 100 150 200 250 300 "Td 50 100 150 200 250 300 

m^^ (GeV) m„ (GeV) 

Figure 3: (Left) Diphoton differential cross section as a function of the photon pair invariant 
mass nijj from data (points) and from theory (solid line) for the photon pseudorapidity range 
\t]\ < 2.5. (Right) The difference between the measured and theoretically predicted diphoton 
cross sections, divided by the theory prediction, as a function of nijj. In both plots, the inner 
and outer error bars on each point show the statistical and total experimental uncertainties. The 
4% uncertainty on the integrated luminosity is not included in the error bars. The dotted line 
and shaded region represent the systematic uncertainties on the theoretical prediction from the 
theoretical scales and the PDFs, respectively. 
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Figure 4: (Left) Diphoton differential cross section as a function of the photon pair invariant 
mass m-y-y from data (points) and from theory (solid line) for the photon pseudorapidity range 
\t]\ < 1.44. (Right) The difference between the measured and theoretically predicted diphoton 
cross sections, divided by the theory prediction, as a function of niyj. In both plots, the inner 
and outer error bars on each point show the statistical and total experimental uncertainties. The 
4% uncertainty on the integrated luminosity is not included in the error bars. The dotted line 
and shaded region represent the systematic uncertainties on the theoretical prediction from the 
theoretical scales and the PDFs, respectively. 
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diphoton cross sections, divided by the theory prediction, as a function of Pt,77- In both plots, 
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tainties. The 4% uncertainty on the integrated luminosity is not included in the error bars. The 
dotted line and shaded region represent the systematic uncertainties on the theoretical predic- 
tion from the theoretical scales and the PDFs, respectively. 
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Figure 6: (Left) Diphoton differential cross section as a function of the photon pair transverse 
momentum pT,yy from data (points) and from theory (solid line) for the photon pseudorapid- 
ity range |/^| < 1.44. (Right) The difference between the measured and theoretically predicted 
diphoton cross sections, divided by the theory prediction, as a function of Pt,77- In both plots, 
the inner and outer error bars on each point show the statistical and total experimental uncer- 
tainties. The 4% uncertainty on the integrated luminosity is not included in the error bars. The 
dotted line and shaded region represent the systematic uncertainties on the theoretical predic- 
tion from the theoretical scales and the PDFs, respectively. 
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Figure 7: (Left) Diphoton differential cross section as a function of the azimuthal angle be- 
tween the two photons, Acp^y, from data (points) and from theory (solid line) for the photon 
pseudorapidity range \rj\ < 2.5. (Right) The difference between the measured and theoretically 
predicted diphoton cross sections, divided by the theory prediction, as a function of Af^^. In 
both plots, the inner and outer error bars on each point show the statistical and total experi- 
mental uncertainties. The 4% uncertainty on the integrated luminosity is not included in the 
error bars. The dotted line and shaded region represent the systematic uncertainties on the 
theoretical prediction from the theoretical scales and the PDFs, respectively. 
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Figure 8: (Left) Diphoton differential cross section as a function of the azimuthal angle be- 
tween the two photons, Acpjy, from data (points) and from theory (solid line) for the photon 
pseudorapidity range |//| < 1.44. (Right) The difference between the measured and theoreti- 
cally predicted diphoton cross sections, divided by the theory prediction, as a function of Acpyj. 
In both plots, the inner and outer error bars on each point show the statistical and total exper- 
imental uncertainties. The 4% uncertainty on the integrated luminosity is not included in the 
error bars. The dotted line and shaded region represent the systematic uncertainties on the 
theoretical prediction from the theoretical scales and the PDFs, respectively. 
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Figure 9: (Left) Diphoton differential cross section as a function of |cos0*| from data (points) 
and from theory (solid line) for the photon pseudorapidity range |?/| < 2.5. (Right) The dif- 
ference between the measured and theoretically predicted diphoton cross sections, divided by 
the theory prediction, as a function of |cos0*|. In both plots, the inner and outer error bars 
on each point show the statistical and total experimental uncertainties. The 4% uncertainty on 
the integrated luminosity is not included in the error bars. The dotted line and shaded region 
represent the systematic uncertainties on the theoretical prediction from the theoretical scales 
and the PDFs, respectively. 
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Figure 10: (Left) Diphoton differential cross section as a function of |cos0*| from data (points) 
and from theory (solid line) for the photon pseudorapidity range |?/| < 1.44. (Right) The dif- 
ference between the measured and theoretically predicted diphoton cross sections, divided by 
the theory prediction, as a function of |cos0*|. In both plots, the inner and outer error bars 
on each point show the statistical and total experimental uncertainties. The 4% uncertainty on 
the integrated luminosity is not included in the error bars. The dotted line and shaded region 
represent the systematic uncertainties on the theoretical prediction from the theoretical scales 
and the PDFs, respectively. 
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Table 3: Measured diphoton differential cross section as a function of pT,jy for the two photon 
pseudorapidity ranges, with statistical (stat.) and systematic (syst.) uncertainties. 



dcr/dpT,^^ [pb/GeV] 



Pt, 


„[GeV] 


\t]\ < 1.44 




\v\ 


<2.5 










Stat. 


syst. 




Stat. 


syst. 


0- 


-4 


0.93 


±0.13 


+0.04 


-0.05 


1.94 


±0.32 


+0.12 


-0.13 


4- 


-6 


1.20 


±0.42 


+0.10 


-0.09 


3.80 


±0.88 


+0.27 


-0.29 


6- 


-8 


1.68 


±0.45 


+0.12 


-0.12 


2.66 


±0.87 


+0.27 


-0.24 


8- 


-12 


1.24 


±0.22 


+0.08 


-0.08 


2.21 


±0.45 


+0.26 


-0.22 


12- 


-18 


0.85 


±0.14 


+0.06 


-0.06 


1.61 


±0.28 


+0.15 


-0.15 


18- 


-30 


0.320 


±0.058 


+0.026 


-0.022 


0.63 


±0.12 


+0.09 


-0.08 


30- 


-40 


0.262 


±0.055 


+0.019 


-0.017 


0.57 


±0.10 


+0.05 


-0.04 


40- 


-50 


0.234 


±0.049 


+0.020 


-0.019 


0.507 


±0.093 


+0.040 


-0.036 


50- 


-80 


0.077 


±0.017 


+0.007 


-0.007 


0.153 


±0.030 


+0.016 


-0.016 


80- 


-180 


0.0084 


±0.0026 


+0.0006 


-0.0005 


0.0150 


±0.0036 


+0.0010 


-0.0009 



Table 4: Measured diphoton differential cross section as a function of Acpjy for the two photon 
pseudorapidity ranges, with statistical (stat.) and systematic (syst.) uncertainties. 



da/ dA^j^ [pb] 





\tj\ < 1.44 


\ri\< 2.5 






stat. 


syst. 




stat. 


syst. 


0-0.271 


1.87 


±0.53 


+0.13 


-0.13 


4.65 


±0.89 


+0.29 


-0.30 


0.2/1-0.471 


1.77 


±0.55 


+0.15 


-0.14 


5.5 


±1.1 


+0.5 


-0.4 


0.471-0.671 


3.09 


±0.72 


+0.31 


-0.29 


5.5 


±1.3 


+0.6 


-0.5 


0.67r-0.87r 


7.2 


±1.1 


+0.5 


-0.4 


16.1 


±2.1 


+1.4 


-1.2 


0.871-0.88 TT 


20.8 


±2.6 


+1.0 


-1.0 


36.7 


±5.3 


+3.4 


-3.0 


0.8871-0.9271 


29.8 


±5.1 


+1.7 


-1.5 


67 


±11 


+5 


-5 


0.9271-0.9571 


36.2 


±8.1 


+5.1 


-4.7 


66 


±15 


+9 


-8 


0.957r-0.987r 


58.8 


±8.8 


+4.2 


-3.8 


103 


±17 


+12 


-11 


0.9871-71 


68 


±11 


+4 


-4 


141 


±23 


+12 


-11 
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9 Summary 



Table 5: Measured diphoton differential cross section as a function of |cos 6* \ for the two photon 
pseudorapidity ranges, with statistical (stat.) and systematic (syst.) uncertainties. 



do-/d\cos6*\ [pb] 



|cosr| 


< 1.44 


<2.5 






Stat. 


syst. 




Stat. 


syst. 


0-0.2 


52.6 


±5.2 


±3.1 -3.2 


87.3 


±9.0 ± 


9.1 -7.9 


0.2-0.4 


38.4 


±4.9 


±3.0 -3.0 


67.0 


±8.2 ± 


6.6 -6.0 


0.4-0.6 


34.8 


±4.6 


±2.7 -2.5 


66.0 


±7.5 ± 


5.9 -5.3 


0.6-0.8 


25.6 


±3.7 


±1.6 -1.5 


66.7 


±7.7 ± 


6.1 -5.3 


0.8-1 


6.4 


±1.4 


±0.3 -0.4 


30.8 


±7.9 ± 


5.9 -4.7 



by the CMS detector in 2010, corresponding to an integrated luminosity of 36 pb^ . The dif- 
ferential cross sections have been measured as functions of the diphoton invariant mass, the 
diphoton transverse momentum, the difference between the two photon azimuthal angles, and 
I cos 9*\. The background from hadron decay products was estimated with a statistical method 
based on an electromagnetic energy isolation variable X. The signal and background distribu- 
tions for I were entirely extracted from data, resulting in systematic uncertainties of approxi- 
mately 10% on the measured diphoton yields. 

The measurements have been compared to a theoretical prediction performed at next-to-leading- 
order accuracy using the state-of-the-art fixed-order computations [1, 2J. Whereas there is an 
overall agreement between theory and data for the diphoton mass spectrum, the theory un- 
derestimates the cross section in regions of the phase space where the two photons have an 
azimuthal angle difference Acp < 2.8. 
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